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ABSTRACT 

A companion paper presents a worked model for evolution through inflation to initial 
conditions for an isocurvature model for structure formation. It is shown here that the 
model is consistent with the available observational constraints that can be applied 
without the help of numerical simulations. The model gives an acceptable fit to the 
second moments of the angular fluctuations in the thermal background radiation and the 
second through fourth moments of the measured large-scale fluctuations in galaxy counts, 
within the possibly significant uncertainties in these measurements. The cluster mass 
function requires a rather low but observationally acceptable mass density, 0.1 <^ 0.2 
in a cosmologically flat universe. Galaxies would be assembled earlier in this model 
than in the adiabatic version, an arguably good thing. Aspects of the predicted non- 
Gaussian character of the anisotropy of the thermal background radiation in this model 
are discussed. 

Subject headings: cosmology: theory — cosmology: large-scale structure of universe — 
galaxies: formation 
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1. Introduction 

An accompanying paper (Peebles 1998a; hereafter 
Paper I) presents a worked example of the evolution of 
a cosmological model through inflation to initial con- 
ditions for an isocurvature (ICDM) model for struc- 
ture formation in a universe that now is dominated 
by cold dark matter. Here I show that the model can 
be adjusted to flt main observational constraints. 

As in Paper I, I attempt to keep the discussion sim- 
ple and definite by adopting a specific set of model 
parameters chosen to give a reasonable approxima- 
tion to the observations. More detailed parameter 
studies that seek to minimize measures of fit to 
the full suite of constraints would be interesting but 
perhaps are not yet a pressing need because many im- 
portant observational constraints still are preliminary 
and may harbor systematic errors. 

The adiabatic cold dark matter (ACDM) model 
for structure formation has been subject to search- 
ing tests from numerical simulations (eg. Governato 
et al. 1998; Springel et al. 1998; and references 
therein). I hope the simpler observational tests pre- 
sented here show that the considerable effort needed 
for a meaningful application of numerical simulations 
of the ICDM model would be worthwhile. 

The model parameters are listed in §2. Second 
moments of the angular distribution of the thermal 
background radiation (the CBR) and the large-scale 
space distribution of galaxies are presented in §3. In 
the ICDM model the primeval CDM mass distribu- 
tion is proportional to the square of a random Gaus- 
sian process with zero mean. In §4 I discuss the na- 
ture of the large-scale non-Gaussian fluctuations in 
the mass distribution and compare them to third and 
fourth moments of galaxy counts. The mass function 
of rich clusters of galaxies is discussed in §5. Because 
the distribution of mass fluctuations is broader than a 
Gaussian with the same standard deviation, rare mass 
concentrations form earlier than in an ACDM model. 
The ICDM model thus requires a lower mean mass 
density for given normalization of the power spec- 
trum, and the cluster mass function changes signif- 
icantly less rapidly with redshift than in the ACDM 
model. In §6 I present the scaling relation between the 
epochs of assembly of the dark matter concentrations 
in galaxies and in rich clusters of galaxies. The rela- 
tively early assembly of protogalaxies in the ICDM 
model is arguably attractive. Finally, §7 presents 
some considerations of the higher moments of the 



angular fluctuations of the CBR. As an example I 
compute the third moments of the quadrupole and 
octupole components of the CBR anisotropy. Con- 
cluding remarks are presented in §8. 

2. Model Parameters 

The cosmological parameters are the same as in 
Paper I, 

f7 = 0.2, f7B;$0.05, A = 0.8, 

To = 2.73 K, Ho = 70 km s"^ Mpc"^ (1) 

The density parameter in baryons is U,^ , and the den- 
sity parameter in CDM is — VI-q . The model is cos- 
mologically flat; A is the fractional contribution to 
the square of the expansion rate by a term in the 
stress-energy tensor that is (or acts like) a cosmolog- 
ical constant A. 

The primeval entropy per baryon is a fixed univer- 
sal value, to agree with the standard model for the 
origin of the light elements. There are no spacetime 
curvature fluctuations at high redshift: the net mass 
density is homogeneous. Homogeneity is broken by 
the irregular primeval distribution of the CDM, which 
is assumed to be a massive scalar field (or its decay 
remnants) squeezed by infiation from its ground level 
to a classical field 0(x). The resulting mass distri- 
bution after infiation ends and before the field starts 
oscillating (when the Hubble length increases to the 
Compton wavelength of the CDM) is 



p(x) = m20(x)V2, 



(2) 



where (/'(x) is a random Gaussian process with zero 
mean and power spectrum 



-2.4, 



(3) 



on scales of interest. The field autocorrelation func- 
tion is 



^P,e-- « (4) 

for the values of of interest here (eq. [||). As 
discussed in §4, the mass autocorrelation function is 



^ = 2 



(0(xi)(/«(xi))2 -(6+2m^) 



A2p 



(5) 
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and the mass fluctuation power spectrum is 

= / d^x^{x)e'^ '' oc (6) 



where 



= 3 + 2mdi = — 1. 



(7) 



The numerical value fits the second moments of the 
CBR and large-scale galaxy distributions (§3). The 
simplest inflation model for isocurvature initial condi- 
tions gives niff, = — 3; the model in Paper I is arranged 
to produce the wanted "tilt" in equations (^) and <^). 
The power spectrum is normalized to 



P{k) 6300/1"'^ Mpc" at k = O.lh Mpc" 



(8) 



where the adopted Hubble parameter is h — 0.7 
(eq. 0). The rms fluctuation in the mass in a sphere 
of radius 8h~^ Mpc is as = 0.9. By this traditional 
measure the model is biased, but as discussed next the 
normalization fits one of the better measurements of 
the second moment of the large-scale galaxy distribu- 
tion. 

Finally, the ionization history is computed in the 
standard way, taking account of the slowing of recom- 
bination by the Lyman-a recombination radiation, 
and under the assumption that there is no significant 
source of ionizing radiation other than the CBR and 
the recombination. 

3. Second Moments of the Distributions of 
Mass and Radiation 

The spherical harmonic expansion of the angular 
distribution of the CBR is 

m<^)=.^arrr(^,'/'), (9) 

and the second moments may be written as 
'1(21 + 1)'^^/^ 



Ti 



4tt 



(10) 



In the approximation of the sum over I as an integral, 
^ / ~ J dl /I, the variance of the CBR temperature 
per logarithmic interval of / is (T;)^]^ 



^ There are good historical reasons, dating from the introduction 
of the ACDM model, for writing 2l{l+l) in place of 1(21+1), but 
since I am considering the ICDM case the use of the convention 
in equation ( |lO[ ) that takes account of the 21 + 1 values of m 
for given / may be reasonable. 



The ICDM model results in Figure 1 use the pa- 
rameters in §2, and show the effect of adjusting the 
baryon density parameter from = 0.05 to 0.01. 
The measured T; are from the compilation of Ratra 
(1998). The model is low at I ~ 100, and perhaps also 
at Z = 2, depending on the correction for the Galaxy. 
Gawiser and Silk (1998) present more detailed com- 
parisons to the measurements. I conclude that since 
these difficult measurements may contain undetected 
systematic errors the model fit is about as good as 
might be expected. 

The power spectrum of the spatial distribution of 
mass in the ICDM model is shown in Figure 2. The 
data are from the IRAS PSC-z (point source cata- 
log) redshift survey of Saunders et al. (1998). This is 
the real space spectrum after correction for peculiar 
velocity distortion represented by the density-bias pa- 
rameter [3 = 0.6. There are good measurements of the 
spectrum of the galaxy distribution on smaller scales, 
k > O.lh Mpc~^, but this approaches the nonlinear 
sector, and it seems appropriate to postpone discus- 
sion of structure on relatively small scales until we 
have more detailed explorations of nonlinear evolu- 
tion from the non-Gaussian initial conditions of this 
model. Since the PSC-z catalog is deep, with good 
sky coverage, it promises to be an excellent probe of 
the large-scale galaxy distribution. Again, the model 
fit seems to be as good as might be expected. 

4. Higher Moments of the Galaxy Counts 

In the ICDM model the primeval CDM mass dis- 
tribution is the square of a Gaussian with zero mean 
value. I consider first some statistical properties of 
the distribution and then compare the moments to 
galaxy counts. 

4.1. Statistical Character of the Mass Fluc- 
tuations 

The mass density is p(x) oc 0(x)^, where {(p) — 
and the autocorrelation function of the Gaussian field 
is 



(0(xi)0(x2)) oc a;i2 



(mp + 3)/2 = 0.6, (11) 



for the power spectrum index nip in equation (0). 
The fluctuation spectrum is suppressed on large scales 
where radiation pressure cannot prevent the tendency 
of the isocurvature model to remain homogeneous, 
and nonlinear evolution distorts the primeval spec- 
trum on small scales. Here I consider intermediate 
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scales, k ^ 0.1 /i Mpc~^, where the shape of the mass 
fluctuation spectrum is close to the primeval form rep- 
resented by equation (pi]). 

The mass distribution smoothed through a window 
of volume V, with the mean subtracted (and units 
chosen so the field mass is unity), is 

Sps = I ^(0(x)^ - (0^)), (12) 

and the n}^ central moment is 

/rl"-V 
— (n(0(xi)2-(02^)). (13) 

The integrand is the expectation value of a sum of 
products of field values at the positions of the vari- 
ables of integration x^. Since 0(x) is a Gaussian 
process the expectation value of each term is the 
sum of products of two-point correlation functions 
(0(xi)(/)(xj)) for all ways of pairing all positions 
All cases with i = j are eliminated by the term {(f>^), 
leaving the terms in which every position appears in a 
two-point function with a different position. That is, 
the integral is over a sum of terms each of which is a 
product of n factors Xij'^ with i ^ j. For < e < 1.5, 
as in equation (0), the integrals converge at small 
Xij and scale with the window size r (the radius, if a 
spherical window) as (((Jps)") oc r^'". In particular, 
the variance scales as {{Sps)'^) oc r^^^. The moments 
of the density contrast S — 5ps/ {p) thus satisfy 

{{6/ar) = F„, {5^'^ (14) 

where Fn is independent of the window size. That is, 
given the window shape the probability distribution 
function of the ratio 5 /a oi the density contrast to its 
standard deviation has a universal form, independent 
of the window size. 

I have not been able to find an analytic expres- 
sion for the distribution function of 5/a for given 
e and window shape, but can report numerical val- 
ues for low moments. The exercise of contracting the 
products of pairs of the (/)(xi) in two-point functions 
yields the reduced spatial two-point correlation func- 
tion (eqs. |], 0), 

6 = {5i52) = 2{xjxi2f\ (15) 

where coherence length, the reduced three- 

point function, 

6 = {5i525:i) - 8r37(xl2X23a;3l)^ (16) 



and the reduced four-point function, 
-{5ih){525z) 

= l&rf[{xi2X2iX^AXAiy + {xi2X2iXA3Xil)~'' 

+ (a;i32;32a;24a;4i)"^]- (17) 
The skewness of the smoothed mass contrast satisfies 

i?3H(V.)^) = 23/^«^f^ggp. (18) 

\^12 I 

The angular brackets on the right-hand side mean the 
integral over the window and divided by the window 
volume, a convenient form for Monte Carlo integra- 
tion. The excess kurtosis satisfies 

D, = {{5/af)--i{{5/a)Y 

,^ ((3:123:233:34x41)"^) , . 

{x^ir ■ ^ ^ 

Numerical values for e — 0.6 and a square (top hat) 
spherical window are 

D3 - (((5/a)3) = 2.46, 

D, = {{S/a)^)-3{{S/a)Y = 9.87. (20) 

4.2. Higher Moments of Deep GalzLxy Counts 

Moments of galaxy counts in catalogs of angular 
positions have been studied to test the theory of the 
onset of the gravitational growth of nonlinear cluster- 
ing out of an initially Gaussian mass distribution. In 
the ICDM model there are initial intrinsic higher or- 
der moments as well as what is generated by the onset 
of nonlinear clustering. This analysis takes account 
of the intrinsic part only. The notation follows Pee- 
bles (1980), with some adjustments to modern con- 
ventions. 

It will be assumed that a galaxy at distance r is 
in the angular catalog with probability proportional 
to a single selection function ilj{r), and that the se- 
lection probabilities along different lines of sight are 
statistically independent. Then the reduced n-point 
angular correlation function is 

w„(l,2,...n) = y^„(l...n) UiJir^dn, 
/ i}{r)r^dr = 1. (21) 

The integral of the reduced spatial function, is 
over the volume element per steradian along each line 
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of sight and weighted by the selection function. The 
angular average of this expression across a field of 
solid angle is 



(22) 



As indicated in the last expression, this is an average 
of the spatial correlation function over n positions 
placed uniformly at random in the angular field and 
with distributions in radial positions given by 



P{< r) 



(23) 



where P is uniformly distributed from P = to P = 
1. Again, this is a convenient form for Monte Carlo 
integration. 

With the correlation functions in equations ( p^ ) 
to (|l^) one sees that average of the three-point angu- 
lar function (eq. |23]) satisfies 



W3_ _ ^3/2 (('^12'-23?-3l) ') 



3/2 



and the average of the fourth moment satisfies 

^^4 _ .„((7'l27'23''34?'4l)"^) 



^4 = ^ = 12 



(24) 



(25) 



These expressions are averages through conical win- 
dows, as in equations ( p^ and (|l^), but here the ra- 
dial distribution is given by equation (p3|). 

The selection function may be modeled in terms of 
a Schechter luminosity function. 



ip{r) oc 



AV^e-V^., (26) 



where the probability a galaxy is in the catalog is the 
probability the energy flux density from the galaxy 
at the observer exceeds the threshold /. This can be 
rewritten as 



ip (X I u 

Ir/D 



l+2a 



e-" du, D, = (i,/47r/)i/^ (27) 



The characteristic depth, i?*, of the catalog does 
not enter the scale-invariant ratios in equations ( p^ ) 
and (H). 

Table 1 shows numerical values of the ratios d^ and 
di (eqs 1 24 and |^^) computed using the selection 
function in equation (E7h with a = — 1.0. At a = —1.5 



the ratios d^ and d4 are about 10% smaller. These 
ratios for narrow cones are smaller than for a spherical 
window (eq. [Q), but the difference is not large. 

The relations between the means Wn of the corre- 
lation functions and the moments of counts in cells is 
discussed in Peebles (1980 §36) and in more general- 
ity by Gaztaiiaga (1994). In the former notation, and 
in terms of the central moments 



fJ-n 



{{N-{N)r), 



(28) 



for counts N in cells, the third and fourth reduced 
moments corrected for shot noise are 

{Nfws = M3-3Ai2 + 2(iV), (29) 
{N)^Wi = ^4 -3^*2 -6^*3 + 11/^2 -6(iV). 

Table 1 shows ratios of the Wi derived from galaxy 
counts in the Edinburgh-Durham Southern Galaxy 
Catalogue (Collins, Nichol, & Lumsdcn 1992) and the 
APM Catalogue (Maddox et al. 1990). The fields are 
square, with solid angle 6x9. In the analysis of the 
onset of the gravitational growth of nonlinear cluster- 
ing out of an initially Gaussian mass distribution one 
is interested in the ratios 



(30) 



I convert the estimates of s„ from the EDSGC catalog 
(Table 1 in Szapudi, Mciksin & Nichol 1996) to the 
dn using 

u,2 = 0.079-°-^'^, (31) 

where 9 is measured in degrees, from a power law fit 
to the two-point correlation function w{9) at 0.1° < 
9 < 2° (Nichol & CoUins 1993). Gaztaiiaga (1998) 
kindly provided the dn derived from his analysis of 
moments of counts in the APM catalog (Gaztaiiaga 
1994). 

Szapudi & Gaztahaga (1998) point out that the 
APM and EDSGC catalogues are independently ob- 
tained from the same photographic plates. They show 
that the moments of counts in the subsample of APM 
in the smaller field of EDSGC are in satisfactory 
agreement with the EDSGC moments. The substan- 
tial difference of the da and d4 at 9 = 2° thus is in the 
sky or plates or more limited size of the EDSGC field. 
Szapudi & Gaztanaga note that since the APM field 
is substantially larger it likely is the more reliable, 
but that can leave room for appreciable uncertainty 
in the moments from APM. 



5 



A simpler measure may be relevant: the two- 
and three-point spatial functions in equations 
and ( p^ ) extrapolated to — xi2 — 123 = x^i in 
linear theory satisfy ^2 = 2, ^3 = 8, and Q = 0.7, 
close to the observed value in the hierarchical form 
for the galaxy three-point function. Thus I conclude 
that the measurements are not inconsistent with the 
model prediction under linear perturbation theory. 

A more serious challenge comes from the correc- 
tion for the nonlinear growth of clustering. In the nu- 
merical analysis of evolution from the non-Gaussian 
initial conditions of a texture model Gaztahaga & 
Mahonen (1996) find the skewness parameter grows 
from the initial value D3 = 0.7 at cr = ((J^)!/^ = 0.1 to 
£>3 = 4 at cr = 1. Comparable growth from the initial 
value £)3 = 2.5 in the ICDM model (eq. ||^) could 
make unacceptably large skewness, but that does de- 
pend on the relative effect of the primeval skewness 
and excess kurtosis on the growth of the variance 
and skewness. Potentially powerful probes of possible 
departures from non-Gaussian initial conditons from 
the properties of clusters of galaxies are discussed by 
Chiu, Ostriker, & Strauss (1998), and Robinson, Ga- 
wiser, & Silk (1998). It remains to be seen whether 
the initial conditions of the ICDM model evolve to 
the observed clustering hierarchy of the galaxy distri- 
bution at (T ^ 1 (Scoccimarro et al. 1998; Fosalba & 
Gaztahaga 1998). 

5. The Cluster Mass Function 

It is usually agreed that the Press-Schechter (1974) 
approximation offers a useful way to estimate the 
mass fraction in rare concentrations such as rich clus- 
ters of galaxies when the primeval density fluctuations 
are Gaussian. Here I apply a variant of the Press- 
Schechter method based on a numerical determina- 
tion of the mass fraction in rare peaks in the ICDM 
model. 

The numerical determination starts with a realiza- 
tion of a random Gaussian process with power law 
power spectrum in a 256'^ cubic grid of positions and 
wavenumbers. For the ICDM model the realization 
is squared and then averaged on the lattice points 
through a square (top hat) window of radius rg to 
get a smoothed mass distribution. Peaks are defined 
as points on the lattice where the density is larger 
than at the six nearest lattice points. The mass 
fraction in peaks is the fraction of randomly placed 
points that fall within distance of a peak with 



smoothed density contrast 8p larger than a chosen 
multiple of the standard deviation, 5p > va. In the 
non-Gaussian case the Gaussian process has power 
spectrum cx /c^^ *, so the mass fluctuation spec- 
trum is Pp oc fc~^'^. In a comparison Gaussian model 
the field is not squared, and here oc k~^^^ . The nu- 
merical results in Figure 3 are based on the smoothing 
window radius = 4.1, in units where the box width 
is 256. Results at half the value of are quite similar. 
An analytic fit for the ICDM model is 

/(> 81a) = 0.37e~"-6^*/'". (32) 

Following the Press-Schechter method, the CDM 
mass distribution at high redshift is smoothed through 
a spherical window that contains the wanted mass, 
and the smoothed density contrast is extrapolated to 
the epoch of interest in linear perturbation theory. 
The mass fraction in peaks in this distribution, with 
density contrast 

(5p = vop > Sc = 1.68, (33) 

is the approximation to the mass fraction present in 
collapsed concentrations at the chosen epoch. Here 6c 
is the critical contrast for spherical collapse at A = 0. 

To get the present mass fraction in clusters I take 
the minimum cluster mass to be 

mci = 4 X IO^^H^^Mq. (34) 

The survey of Bahcall & Cen (1993) indicates the 
present number density of clusters at least this mas- 
sive is 

nci = (2 ± 1) X 10"^/i"^ Mpc"^ (35) 

The cluster mass function varies roughly as n(> m) cx 
TO^^, so the mass fraction in clusters at to > md is 

/ci = 2TOcinci/(p), (36) 

where (p) is the cosmic mean mass density. These 
numbers give 

0.003 < fciQ < 0.009. (37) 

The mass in equation (|3j) is contained in a sphere of 
comoving radius 

Tci = 7.0h-^n-^/^ Mpc, (38) 

in the original near homogeneous mass distribution. 
The mean square fractional fluctuation in galaxy 
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counts in a sphere of this radius in the present galaxy 
distribution (and ignoring shot noise) is reasonably 
well approximated by the power law model 

Ug^ = 1.82[(4.5±0.5) Mpc//irci]^-^^ 
dg ^ (0.91 ± 0.09)f7" '''". (39) 

To test the method I consider first the Gaussian 
case at redshift z — Q. For Vl — 1 the mass frac- 
tion in equation ( |37| ) with the lower curve in Figure 3 
translates to 3.1 < v < 3.5. The collapse condi- 
tion in equation (|3^) yields the standard deviation 
in the mass contrast in the window that contains md, 
0.47 < CTp < 0.54. The ratio to the rms fluctuation in 
galaxy counts from equation (^ is 

0.5 < ap/ag < 0.7. (40) 

For = 0.2 the same calculation for the Gaussian 
model gives 

0.9 < cfp/ag < 1.4. (41) 

Equations 1^) and (|l|) are close to the usual es- 
timates of the bias needed to produce a reasonable 
cluster mass function at the present epoch (Eke, Cole, 
& Frenk 1996; Cen 1998). That is, this apphcation 
of the Press-Schechter procedure based on the mass 
fractions in peaks in a simulation seems to be reason- 
ably secure for the Gaussian case. 

For the non-Gaussian case the procedure using the 
upper curve in Figure 3 with ft — 1 gives 

0.2 < ap/ag < 0.4. (42) 

For n = 0.2, the value adopted in the ICDM model, 
the method gives 

3.K 1/ < 4.7, (43) 

to fit the mass fraction in clusters, and 

0.6 < ap/ag < 1.0. (44) 

The more extended tail allows mass fluctuations at 
larger values of v = S/a than in the Gaussian case. 
A result is that if = 1 the suppression of the fluc- 
tuations in mass relative to galaxies has to be even 
stronger than in the Gaussian case. For fl = 0.2 equa- 
tion (^) indicates the ICDM model is just within the 
bound from the cluster mass function and the assump- 
tion galaxies trace mass, ap ~ ag. 

The predicted evolution of the cluster mass func- 
tion is much slower in the non-Gaussian model than 



in the Gaussian case because the broader tail of den- 
sity fluctuations makes the mass fraction a less sen- 
sitive function of S/a. For example, for fl = 0.2 the 
growth factor for the density contrast in linear per- 
turbation theory from redshift z = 0.5 to the present 
is D = 1.24. When v in equation (^) is multiplied by 
this factor it reduces the mass fraction /d (from the 
top curve in Fig. 3) by a factor of two, meaning the 
comoving cluster mass function at z = 0.5 is half the 
present value. The evolution is faster if = 1: the 
present central value v = 6.2 translates to v = 9.3 
ai z — 0.5, and the corresponding mass fraction is 
/ci ^ 0.001, about 10 percent of the present value. 
This is not inconsistent with the observations of clus- 
ters at z ~ 0.5, but the strong bias is not attractive 
(eq. 11). 

6. The Epoch of Galaxy Assembly 

In the CDM family of models, galaxies are assem- 
bled as mass concentrations by a scaled version of the 
assembly of the mass present now in rich clusters of 
galaxies. One sees from Figure 2 that in linear pertur- 
bation theory the mass fluctuation power spectrum in 
the ICDM model is close to the primeval power law 
form at fc ^ O.lh Mpc~^. Here the spectrum up to 
the onset of development of nonlinear structure varies 
as 

Pp oc /fc™''D(<)2, (45) 

where D{t) is the solution to the linear equation for 
the evolution of the density contrast 5p/p. The rms 
contrast through a window of comoving radius x thus 
scales as 

Structure formation is triggered by passage of upward 
fluctuations through 5^ ~ 1, meaning the comoving 
length scale on which structure is forming varies with 
time as 

Xni CX 1)2/(3+"^). (47) 

The corresponding physical length varies as 

r„i oc (l + 2)-iD2/(3+m,)^ ^4g^ 

the characteristic masses of newly forming objects is 
TO„i ocD'^/t^+^p), (49) 

and the characteristic velocity dispersion within de- 
veloping structures is 

a„ioc(H-z)i/2i?2/(3+m,)^ (50) 
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These relations may be normalized to the great 
clusters of galaxies, with 

TA = 1.5/i~^ Mpc, cTci = 750 km s"\ (51) 

The line of sight velocity dispersion is an rms mean 
for R > I clusters. In the limiting isothermal sphere 
model these numbers yield the mass in equation (|3^). 
Clusters are still relaxing at the Abell radius ta, and 
the merging rate is not insignificant, but it is thought 
that the internal velocities typically are close to what 
is needed for support against gravity at r ta- In 
the power law model in equation ( ^5|) these quantities 
scaled back in time using equations ( ^8|) to (|5^) char- 
acterize objects in a like state of early development in 
the past. 

With the parameters in equations (Q) and (0) the 
scaling relations applied at expansion factor 1 + z = 7 
give 

Tg = 20 kpc, cTg = 140 km s~^, 

mg = 2 X 10"Mo. (52) 

This scaling calculation ignores the dissipative set- 
tling of the baryons to form the luminous central parts 
of the galaxies. Consistent with this, in the standard 
model for an L ^ galaxy the mass within the ra- 
dius r ~ 20 kpc is dominated by dark matter. 

In the ICDM model the large objects at 1 + z = 7 
with the parameters in equation ( |5^ ) may be com- 
pared to present-day clusters: there is significant sub- 
structure, the mass distribution in the outer parts is 
disordered, and there is a significant rate of merging, 
but the internal motions typically are close to what is 
needed for virial support at radius r ^ rg. The mass 
fraction in these objects at 1 + z = 7 is the same as 
the cluster mass fraction from which they are scaled, 
/g - (0.006 ±0.003)/17 ~ 0.03 ±0.015 (eq. ||l). The 
product of /g with the mean mass density and di- 
vided by TTig gives a characteristic comoving number 
density, 

rig = (0.004 ± 0.002)/i3 Mpc~^ (53) 

The characteristic size, internal velocity, mass, and 
number density of these newly assembled systems are 
roughly typical of present-day L* galaxies. An ob- 
server might be inclined to call them young galaxies, 
assembled at z ~ 6. 

The scaling picture says the matter at 1 + z = 7 
that is not in this generation of proto-galaxies would 



be in smaller clouds between them, maybe positioned 
to become the Lyman-a forest observed at z < 5. 

At expansion factor 1 -|- z = 20 the scaling relations 
give 

r — 1.3 kpc, cr = 40 km s~^, 
m=lxlO'^MQ, (54) 

numbers characteristic of dwarf galaxies. I have to 
postulate that some of these objects merge to con- 
tribute to the mass concentrations near the luminous 
parts of present-day giants. In the model the gener- 
ation oi L ^ Lf, protogalaxies has been assembled at 
close to the present comoving number density by the 
epoch 1 -I- z = 7 (eq. |Q ) , so I must assume the rate 
of merging decreases, perhaps because the dissipative 
settling of the baryons has progressed far enough to 
lower merging cross sections. Thereafter structure 
formation would build the present-day galaxy cluster- 
ing hierarchy, while adding the extended massive ha- 
los of galaxies at r ~ 200 kpc by accretion of relatively 
low mass star clusters and diffuse material. These 
numerous postulates could and should be checked by 
numerical simulations. 

Structure formation in the ACDM model is qual- 
itatively the same — hierarchical growth by gravity 
— but it happens at lower redshift (Kauffmann 1996 
and references therein). The remarkable advances in 
the observations of high redshift young galaxies or 
their precursors may be bringing us close to a test 
of these two picture. If galaxies were assembled as 
mass concentrations at 1 -I- z = 7 then the observa- 
tions at 1 + z ~ A ought to reveal internal velocities 
characteristic of present-day L ^ galaxies. This is 
not inconsistent with the properties of the damped 
Lyman-a absorbers studied by Wolfe & Prochaska 
(1998), though Haehnelt, Steinmentz & Ranch (1998) 
show other interpretations are possible. The expected 
optical appearance of young galaxies at 1 -I- z ~ 4 de- 
pends on how feedback affects the rate of conversion 
of gas to stars, a delicate issue that will require in- 
formed discussion. 

7. Higher Moments of the Thermal Back- 
ground Radiation 

This discussion is limited to general remarks and 
an example, the computation of the third moments of 
the quadrupole and octupole parts of the anisotropy 
of the CBR predicted by the ICDM model. 
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Because the statistical properties of the primeval 
mass distribution in the model are more easily ex- 
pressed in terms of the position correlation functions 
than the power spectrum and higher moments of the 
Fourier coefficients, I use a Greens' function repre- 
sentation of the relation between the primeval fluc- 
tuations in the CDM distribution and the observed 
CBR anisotropy. With periodic boundary conditions 
the primeval CDM mass distribution is 



<5(x) =^5ke' 



k X 



(55) 



and the spherical harmonic components of the ob- 
served CBR angular distribution are (eq. ||]) 



^<5ke,(fc)rr(k). 



(56) 



Here e/(fc) is the Legendre transform of the angular 
distribution of the CBR produced by a single Fourier 
component with wavenumber k and normalized to 
amplitude |(5k| that is independent of k. The spherical 
harmonic is evaluated at the direction of k. On ex- 
pressing i5k in equation (|5^) as an integral over S{x), 
and using 

e*k-x = 47r ^ i^ji {kx)Yi'^ (k)>^"'" (x) , (57) 



one can reduce equation (|5^) to 



d^xdix.)Wiix)Yi"\±) 



(58) 



Here the spherical harmonic is a function of the di- 
rection of the position x, and the weight function is 



Wz(a;) = ^y^ edkei{k)ji{kx) 



(59) 



The low order correlation functions of (5(x) in the 
model are given by equations ( p^ ) to ([l^). Equa- 
tion (|5^) expresses a™ as an integral over 5(x), so it 
is straightforward to write down expressions for the 
moments of the a[" as integrals over the spatial cor- 
relation functions. 

The third moments of the a™ may be expressed in 
the form 



\ II l2 l3 ' 



mi 12 



)(i«ri^ 



))l/2 



N 
D' 
(60) 



The numerator is 

N^&J Sxi(fx2d^xz (siaxsssgi)"' 
Wi, {l)Wi, (2) W^,3 (3)1^;-^ (1)1^™^ {2)Yll^^ (3), (61) 
and the denominator is 

D = {Di,Di,Di,f/\ (62) 

where 



Di = 2 



r 

/ d^xid^X2 Xi 



2e 



Wiil)Wii2)Yr{l)Yr"'{2). (63) 

Magueijo (1995) and Ferreira, Magueijo, & Gorski 
(1998) point out that the ratio of moments of the a™ 
in equation ( |60| ) is a sensible way to normalize this 
measure of the departure from Gaussian fluctuations 
in the CBR. The ratio also very conveniently elimi- 
nates the normalizations of the two- and three-point 
mass correlation functions, as in equation (|2j) for the 
skewness of the distribution of galaxy counts. 

The integrals over angular positions are simplified 
by the assumption that the density fluctuations are 
a stationary random process, meaning the correlation 
functions of the (5(x) can only depend on relative po- 
sitions. Thus in the expression for Di one can average 
over all directions at fixed angular separation of the 
lines of sight (1) and (2) using (Peebles 1973) 



(yr(l)V"(2))=^^Kcos^?i2)/4^. 



(64) 



(It would be even easier to go back to the expression 



for 



') as an integral over the power spectrum. 



but it seems best to compute numerator and denom- 
inator in eq. [ |60| in the same way.) 

Luo (1994) shows one can similarly simplify the 
angle integrals over the products of three spherical 
harmonics in equation (|6ll). Because the Yf^ have 
parity (—1)' it immediately follows that nonzero N 
requires 

+ ^2 + ^3 = even integer. (65) 

For even sums the symmetry of rotation about the 
z-axis for given angular separations of the three di- 
rections Xi, X2, and X3 requires 



TOi + m2 + m3 = 0. 



(66) 



As in quantum mechanics, the symmetry of rotation 
about any other axis requires the triangle rule. 



I/1-/2I </3<^l+^2. 



(67) 
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To see this explicitly one can use the method in Pee- 
bles (1973, eq. [B5]), which follows Edmonds (1957), 
to find the mean of the product of three spherical 
harmonics averaged over orientations for given angu- 
lar separations of the three lines of sight, 

(i^7^(i)y,^^(2)r™^(3)) = zi,,,,,w^f±^^^ (68) 



where 



>^/?(o)>^r(^i2,o)r,;"(ei3,02 

The Wigner 3j symbols. 



mi ,m2 ,m3 



h h h 
mi m2 m3 



(69) 



(70) 



are proportional to Clebsch-Gordan coefficients, which 
vanish unless the conditions in equations ( |66| ) and ( |6^ ) 
are satisfied. 

In the isothermal CDM model the primeval fluctu- 
ations in the CDM are balanced by opposing fluctua- 
tions in the radiation, so the positive skewness in the 
CDM distribution is reflected in a negative skewness 
in ST/T. The effect on the distributions of the a™ 
may be illustrated by the expectation values of third 
moments for Z = 2 and / = 3. Equation (^8|) (with 
tables of the Wigner symbols, e.g. Edmonds 1957 Ta- 
ble 2) indicates the nonzero skewness coefficients for 
I = 2 in equation (60) satisfy 



0,2, 



0,1, 



Si 



0,0,0 



S. 



2,2,2 

1,1,-2 
2,2,2 



(2/35)i/2r2,2.2, 

(l/70)l/2r2,2.2, 

-(2/35)i/2r2^2,2, 

-(9/105)1/2T2^2,2. 



(71) 



These moments are unaffected by a permutation of 
the rui (eq. |^^). A sign change, mi — > —mi for all i, 
produces the complex conjugate of the product of ob- 
served a™ and does not affect the expectation value of 
this product. That is, one compares the expectation 
values to the real parts of the products of a™. A mea- 
surement of the al" from a full sky observation thus 
yields four measures of T2_2.2, one of which is the real 



part of (aj''^ 
^2„2 



i2 , which is the same as the real part of 
re 

I = 2 and I — 3 components are 



a 



At 

2 )^a2- The nonzero third moments involving the 



S'. 



3,3,2 



(2/105)1/2T3,3,: 



^1,-1,0 



3,3,2 

3,-3.0 
■'3,3,2 

0,-1,1 
3,3,2 



r.1,-2,1 
•^3,3,2 

r.2,-3,1 
•^3,3,2 

r,-l,-l,2 
•^3.3,2 



s. 



3,3,2 

1.-3,2 
3,3,2 



-(3/140)1/2^^,3,2, 

(5/84)i/2r3,3,2, 

-(1/210)1/2^3,3,2, 

(i/28)i/2r3,3,2, 

-(5/84)l/2r3,3,2, 

(2/35)i/2r3,3,2, 

-(i/2i)i/2r3,3,2, 

(l/42)l/2r3,3,2. 



(72) 



These are nine measures of 13,3,2. 

A Monte Carlo numerical integration of equations 
( |6l| ) and (^ with equations (^) and ( |68| ) is straight- 
forward (and the number of trial positions, ~ 10*, 
needed for convergence no longer a problem). It is a 
comforting check that the results for r332 from — 2, 



Z2 = 3 and from li — 3 and Z2 = 2 in equation (69) 
agree. I find 



1.02, 



T: 



3,3,2 



-0.87. 



(73) 



The analysis of third and fourth moments of other 
combinations of low I multipole components is straight- 
forward in principle, tedious to do by hand, perhaps 
easily done by computerized algebra, but beyond the 
scope of this paper. I expect that for significantly 
larger values of I the best approach will be to analyze 
small sections of the sky in rectangular coordinates 
that allow a much simpler representation of the higher 
order moments. 

In pioneering analyses Kogut et al (1996) found 
constraints on some models for non-Gaussian fluctu- 
ations in the CBR angular distribution, and Ferreira, 
Magueijo, & Gorski (1998) found a set of measures 
of third moments of the aj" in the COBE DMR sky 
maps. I refrain from a comparison to equation (73) 
because the abbreviated moment analysis presented 
here has little overlap with the measures obtained by 
Ferreira, Magueijo, & Gorski (1998). 

Because the initial mass distribution is homoge- 
neous in the isocurvature ICDM model there is no 
significant Sachs-Wolfe effect at I <^ 30; the CBR 
anisotropy simply reflects the non-Gaussian fluctua- 
tions in initial composition. At smaller scales radi- 
ation pressure is able to rearrange the net mass dis- 
tribution, producing significant contributions to the 
CBR anisotropy from the Sachs- Wolfe effect and the 
motion of the baryons. This rearrangement shifts 
phases of Fourier components; thus a prominent hole 
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in the initial distribution of baryons and radiation be- 
comes a ripple. I suspect this significantly suppresses 
the ICDM prediction of non-Gaussian fluctuations in 
the GBR at I ^ 100. 

8. Discussion 

8.1. Is the ICDM Model Attractive from a 
Theoretical Point of View? 

Paper I compares the fields, parameters, and func- 
tional forms of the potential energy in the ICDM 
model and other physical theories for the seeds of 
structure from inflation. Here I consider some broader 
issues. 

Our experience in particle physics might lead us 
to suspect that the laws of physics relevant to the 
early universe will be found to be elegant and simple, 
albeit in some deeply subtle way, and that once we 
understand the physics we will see that the universe is 
an expression of the physics. This world view informs 
many studies of inflation. The night thought of a 
physical scientist might be that Newtonian mechanics 
is expressed in fully developed turbulence, but it is 
not likely we would know much about turbulence if 
we had not seen it. A knowledge of the physics of 
the early universe might not be of much use if its 
expression were complex. 

Two themes could be accepted in either world view 
(as well as by those, perhaps the majority, with more 
moderate opinions). First, the construction of a spe- 
ciflc internally consistent example of how evolution 
from very high rcdshift could have led to the present 
state of the universe is a valuable demonstration of 
consistency of the set of ideas on which it is based. 
We have examples from the adiabatic CDM family 
of models. I have argued for yet another, an isocur- 
vature GDM model. Second, the models and their 
parameters will be reconsidered with each significant 
advance of knowledge of the physics and astronomy, 
a process that will lead us to abandon some models, 
adjust others, and maybe introduce new ones. Per- 
haps this process will back us into that narrow corner 
of model and parameter space that is a useful ap- 
proximation to what really happened. We may have 
a modest example in the fact that this latest version 
of the isocurvature model (earlier steps of which may 
be traced back through Peebles 1997a) has the same 
dynamical actors as the ACDM model, though it re- 
mains to be seen whether this is a lasting situation. 

A mature physical theory must be falsifyable; there 



is good reason for our conditioned dislike of theories 
that can be adjusted to fit whatever is measured. On 
the other hand, if the evolution of the early universe 
were moderately complex wc likely would need a flex- 
ible model to flt it. The isocurvature IGDM model 
in Paper I assumes power law inflation because that 
makes it easy to select the fields and their poten- 
tial energy functions to produce a power law CDM 
fluctuation spectrum Pp oc k^p over a wide range of 
scales. But the evidence may lead us to another func- 
tional form. In the power law IGDM models shown 
as the solid lines in Figure 1 the GBR anisotropy Ti 
at I ~ 100 may be too high (Netterfield et al. 1997). 
That could be remedied by taking the power law in- 
dex nip to be closer to zero, but that would make T; 
unacceptably small at / <^ 10. The dotted line shows 
one way out: change the power spectrum to 

Pp oc {k^/kf + {k,/kf/\ (74) 

with 

k^ = O.Olh Mpc-\ (75) 

and the other parameters the same as for the middle 
solid line. There has to be another bend to P oc fc*, 
s < —3, at fc 1 pc^^ (Paper I). These bends are 
quite inelegant, unless Nature has chosen them. 

The flexibility of the IGDM model is limited. For 
example, it is difficult to lower the spectrum at I ~ 
100 without significantly lowering the peak at I ^ 300. 
The advances in observational constraints from work 
in progress will show whether the IGDM model is a 
useful approximation. If the observations to come in 
the next decade are fitted in all detail by one of the 
simple structure formation models now under discus- 
sion it will compel acceptance. If improving obser- 
vations require increasingly baroque models it may 
mean we have missed the correct elegant picture, or 
that the cn'olution of the universe does not agree with 
our standards of elegance. 

In my reading of the first of the world views men- 
tioned above the ICDM model is quite inelegant be- 
cause it was constructed ad hoc to fit the observations 
and it is flexible enough to be capable of adjustment 
to fit some substantial changes in the observational 
situation. In the second world view, a model that fits 
significant observational constraints within a sensible 
reading of the physics may not be all bad. 
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8.2. Is the Model Attractive from a Phe- 
nomenological Point of View? 

The ICDM model has some possibly significant 
successes and problems. Both will have to be re- 
considered with each advance of the observations and 
their interpretation, of course. I hope is is not entirely 
self-serving to note that a model that is close to real- 
ity may encounter apparent problems as we sort out 
the ambiguities in the evidence. Following arc some 
considerations. 

i) The model parameters that fit the CBR angular 
fluctuation spectrum Ti in Figure 1 fit the second 
moment P{k) of the large-scale galaxy distribution 
in Figure 2, a significant success. The low measured 
value of at Z ~ 100 may require adjustment of the 
model, or perhaps will prove be in some part a sys- 
tematic error in exceedingly difficult measurements. 

ii) The skewness and excess kurtosis of galaxy counts 
are not far from that of the model, a not insignificant 
result. The major open issue is the correction for 
nonlinear evolution: does the non-Gaussian primeval 
mass distribution of the model evolve into the galaxy 
clustering hierarchy? An example of the predicted 
non-Gaussian higher moments of the multipole ex- 
pansion components a™ of the angular distribution 
of the CBR is presented in §7. A comparison to the 
measurements remains to be done. 

iii) The cluster mass function agrees with a Press- 
Schechter approximation under the assumption that 
galaxies trace mass at $7 = 0.2. Chiu, Ostriker, & 
Strauss (1998) point out that this probes the nature 
of the primeval mass fluctuations, because the cluster 
mass function depends on the tail of the distribution 
and the rms galaxy peculiar velocity field on the stan- 
dard deviation. The density parameter in the ICDM 
model discussed here, ^ = 0.2, agrees with the low pe- 
culiar velocities indicated by many analyses (e.g. Pee- 
bles 1986; Bahcall, Lubin, & Dorman 1995; Peebles 
1997b; Wilhck & Strauss 1998). The ICDM model 
thus seems to pass the Chiu et al. test. If further 
work showed that the mass fraction /d in clusters is 
not near the upper end of the range in equation (|37|), 
and the Press-Schechter method is a good approxima- 
tion, it would require a lower value of fi, increasing 
Ti and tending to spoil the general consistency with 
the measured Ti and Pp{k). 

iv) The model predicts relatively early galaxy assem- 
bly; it is an open issue whether this is a success or 
problem. The model may be considered a success 



from a theoretical point of view because I arrived at 
it by a search for galaxy formation at high redshift, 
when the mean mass density would have been consid- 
erably closer to the relatively high density characteris- 
tic of the luminous parts of normal galaxies. The line 
of thought originated in Partridge & Peebles (1967); 
a recent version is in Peebles (1998b). The scaling 
arguments in §6 suggest the ICDM model has some 
attractive features as a model for galaxy formation. 
A more detailed examination by numerical simulation 
remains to be done. 

My conclusion, from the second of the world views 
presented at the beginning of this section, is that the 
ICDM model is attractive because it fits a significant 
set of observational constraints within what appears 
to be an acceptable physical model. 

I am grateful to Bill Ballinger and Will Saunders 
for the PSC-z data in Figure 2; to Neta Bahcall for 
guidance on the cluster data; to Renyue Cen for his 
work on numerical simulations that led to the analysis 
of the cluster mass function in §5; to Avery Meiksin 
for help in interpreting the EDSGC data in Table 1; 
to Enrique Gaztanaga for the APM data and stimu- 
lating discussions of the interpretation of the higher 
moments of the galaxy counts; and to Bharat Ratra 
for the data compilation in Figure 1 and to him and 
Lyman Page for help in its interpretation. This work 
was supported in part by the National Science Foun- 
dation. 
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TABLE 1 
MOMENTS OF GALAXY COUNTS 
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Fig. 1. — Power spectrum of the GBR. The data 
points are from the compilation by Ratra (1998). 



This 2-column preprint was prepared with the AAS lATJiiX 
macros v4.0. 
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The ICDM model assumes the parameters in equa- 
tions (|l|) and (0). The density parameter in baryons 
is Hb — 0.05, 0.03, and 0.01 from the top to bottom 
sohd curve. The dotted curve is discussed in §8. 




k/h (radians Mpc"') 



Fig. 2. — Mass fluctuation spectra extrapolated to the 
present in linear perturbation theory for the model in 
Fig. 1. The data are from the PSC-z collaboration 
(Saunders et al. 1998). The density parameter in 
baryons is flB = 0.05, 0.03, and 0.01 from top to 
bottom at small wavenumber. 




Peak Height i'= 6/a 

Fig. 3. — Mass fractions in peaks at density contrast 
greater than v standard deviations for a Gaussian 
model (lower line) and the non-Gaussian ICDM model 
(upper line). The models have mass power spec- 
tra P (X fc-i-^. The mass distribution is smoothed 
through a spherical window, and / is the mass frac- 
tion within windows centered on peaks of the mass 
distribution smoothed through the window. 
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